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Abstract. The question of interest in the present study is: "Given a body subject 
to mechanical loads, how to define the initial geometry so that the deformed one 
matches precisely a prescribed shape?" This question is particularly relevant in forming 
processes where the tolerated mismatch between the deformed and desired geometries 
may be lower than a Micron. The method proposed here uses as a first "guess" to the 
required initial geometry the desired one, then it updates iteratively the locations of a 
set of boundary points so that their locations in the deformed configuration come closer 
and closer to the desired ones. The scheme is shown to converge unconditionally for 
small deformations in the sense that arbitrarily small mismatch between the deformed 
and desired shape can be achieved. Moreover, since it is based entirely on geometric 
considerations, the convergence should not be affected by the nature of material, i.e. 
it is independent of the constitutive law. The success of the method is illustrated by 
considering an example. 
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1. Introduction 

The manufacturing of high-precision glass or metal products is a challenge due to the 
complexity of the phenomena involved during the forming process and the difficulty to 
control in practice all the relevant parameters. A question engineers often have to face is: 
"given a set of loads applied to a solid body, what should be the initial geometry so that 
the end geometry (after deformation) matches the prescribed one?" The development 
of efficient and accurate numerical techniques such as the Finite Element method 
allowed significant progress towards the solution to this inverse problem. A popular 
approach consists in parameterizing somehow the contour of the work-piece, defining 
an objective function related to the mismatch between the deformed and desired work- 
piece geometries and applying the extensive theory of optimization in order to reduce 
the objective function. This approach has proven its applicability for the problem of tool 
design in the metal forging process (see P3 El and reference therein). Many authors 
favour "gradient-based" optimization techniques and the condition to the success of 
these is the accurate and efficient evaluation of the sensitivities. The sensitivity analysis 
may be performed using the Finite Difference method, the direct differentiation method, 
[21 IS], or the adjoint-state method, [HE]. In many practical cases, the pragmatic user 
of a commercial code for the computation of the deformations is reduced to the use of 
the Finite Difference method to compute the sensitivities and this approach is known 
for its lack of accuracy and efficiency, [Tj. 

An alternative approach which overcome this issue was first introduced by Park et al 
, [3] . It consists in starting from the desired work-piece geometry and back-tracking the 
deformation path. The idea was applied with success by various authors to preform or 
die design in metal forging, [Sid- A similar concept was exploited in jH] in combination 
with an iterative scheme to update the locations of a set boundary points in order 
to identify the required initial geometry. This method proved its applicability and its 
efficiency in the context of high-precision glass forming and was easily implemented 
in a commercial code. In spite of the non-linear viscoelastic behaviour of the glass, 
only few iterations were required to reduce the mismatch between the deformed and 
desired work-piece geometries below a Micron. A formal analysis explained the good 
convergence properties of this iterative scheme and revealed a weakness. It is, in theory, 
only applicable to deformations dominated by shear or volumetric strains. 

The aim of the present work is to present a modified scheme which circumvent this 
restriction. The main idea remains essentially the same but a corrective term is added 
for the update of the boundary points locations. The next two sections are devoted to 
the description of the scheme and the analysis of its convergence properties. The scheme 
is next applied to a case where the prescribed displacement field combines shear and 
volumetric strains in equal magnitude in order to assess its performance. 
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Figure 1. Sketch of the general forward mechanical problem. 
2. Description of the iterative scheme 

We restrict the analysis to a mechanical problem in a two-dimensional domain Q such 
that the displacement field U is a two-dimensional vector whose components depend only 
on the coordinates (x, y) of the Cartesian coordinate system associated with the body 
shown on Figure H The body is subject to mechanical loads. On T u , displacements are 
prescribed and traction is imposed on T a . Clearly, the forward problem corresponds to 
calculating the deformed shape given the initial one and the boundary conditions while 
the inverse problem we wish to solve is: "given a desired shape, find the initial one such 
that the deformed work-piece geometry matches the desired one". We will assume in 
the following that the solution to the forward problem is obtainable by analytical or, 
more likely, numerical means. Consequently, the exact or approximated displacement 
field can be found for all points in Q. 

In order to describe the algorithm which tackles the inverse problem of identifying 
the required initial geometry, notations are first detailed. A number of points are 
introduced at the boundary. If a discrete method is used to solve the forward problem, 
these points naturally correspond to the boundary nodes of the mesh. Let Mf, . . . , Mf 
denote the L boundary nodes of the desired work-piece geometry and M[ m , . . . , M % ™ 
denote the associated L boundary nodes of the required initial geometry (solution to the 
inverse problem). The M™ ,% are found iteratively as the limit of a series of M/, where 
the superscript j denotes the iteration number. At each iteration, N- corresponds to 
the new location of the node M\ in the deformed geometry and U- is the associated 
displacement. The algorithm is best described in pseudo-code notation as follows: 

(i) First iteration : use Mf , i G [1, L] as a first guess to M\ m , i £ [1, L] and calculate 
the residual vector Aj , i G [1, L] whose norm gives a measure of how far the node 
in the deformed geometry is from the desired location: 
for i=l to L { 



6m\ = OM- ; 



(1) 
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ON] = Om) + Ul ; (2) 
Si = OM* - ON] ; (3) 

} 3=2; 

(ii) Following iterations: update the previous estimate to Mf 11 , i G [1, L] and calculate 
the residual vector Aj , i G 
Do { 

for i=l to L { 

OMj = OM j r l + Ar 1 + £| ; (4) 
0%=0M\ + XJ{- (5) 
Al =OM d -0%- (6) 

} 3=3+1;} 

While max(||A^||) > e 

Restated in simple terms, the initial guess for the required initial boundary node 
locations is chosen to be the location of the nodes of the desired geometry. At each 
iteration the residual vector (A;-) whose norm measures how far the deformed geometry 
is from the desired one is evaluated and added to the previous guess of the required 
initial boundary node location. B\ is a corrective term which enlarges the range of 
applicability of the scheme and shall be defined subsequently. If this corrective term 
vanishes, the intuitive scheme proposed in [Hj and illustrated on Figure El is recovered. 
Since this scheme is entirely based on geometric consideration, it does not depend on 
the material behaviour and should therefore be applicable regardless of the constitutive 
law. This was shown in [8 a where the scheme was used for the identification of the 
required initial geometry of a glass piece, a strongly nonlinear material with stress and 
structure relaxation. 



3. Convergence analysis 

In order to assess the convergence properties of the scheme (i.e. will N\ , i G [1,L] 
converge to Mf , i G [1, L\ for increasing j, and if so, how fast?), we perform the same 
analysis as that presented in jH]. Firstly, it is necessary to note that as a consequence 
of the algorithm, 

M = Uf 1 ~U 3 i - B{ . (7) 
This result is proven by simple algebraic manipulations as follows; 
Using eq. (jOJ gives 

K{ - at 1 = oh d -on\- (om* - Oivf 1 ) • 

Substitution of eq. (HJ) into eq. (j3J) and of eq. (0) into the previous equation yields 

M - at 1 = om\~ 2 + Ar 2 + bi 1 + ur 1 - (oMf 1 + ^r 1 + B* + uf) , 
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Figure 2. The iterative algorithm with the corrective term B\ set to zero. 



o- 



but since according to eq. (jH), 0~m\ = OM\ + A{~ 2 + -Bf -1 , eq. (J7J) is recovered 
which proves the result. 

At iteration (j — 1), the displacement of the node M{~ x with coordinates 
{x\~ , Hi' 1 ) can be found by Taylor series expansion around Mf with coordinates (xf, yf). 
Accordingly, 

ujt 1 = u x (xt yf) + (*r - xf) + 9 -^\t (yr - yf) + h. .t. m 

Uyl 1 = U y (xf, yf) + d -^\t (xl 1 - xf) + d -^L\t (yr 1 - yf) + h.o.t. ,(9) 

where {u x {~ 1 ^ U y {~) are the components of the displacement vector U-" 1 and it is 
understood that the spatial derivatives of the displacement are taken at Mf. 
Similarly, at iteration (j), the displacement of the node Mf can be written as follows: 

Ud = U x (xf, yf) + 3 -^\t (4 - xf) + (yl - yf) + h.o.t. , (10) 

Uyl = U y (xf,yf) + (x\ - xf) + ^\f (y\ - yf) + h.o.t. . (11) 

Subtracting eq. (fTUj) to eq. (jHJ) and eq. (fTT|) to eq. (jHJ) gives, 

r.,4 ; - u x { = ^\f (xl 1 - xf) + ^1 (yr 1 - y{) + h.o.t. , (12) 

u v r l - u yi = ^ - 4) + |f I? (yl 1 - vl) + h.o.t. . (13) 

Moreover, taking eqs. (J3J) and J7|) into account and remembering that ^ L = e xx and 
-q£- = e yy , these equations may be rewritten as, 

A i = ... , 'l ( \ J 1 . I ! ' ) - 

dy 



x J , = - <,;! (A,.j 1 + bJ>) ■ ■ — ^ 1 ^r 1 + />',:) - b x \ + h.o.t. , (i n 



corresponds to deformations dominated by volumetric strains (i.e. <C e xx < 1 and 
<C t yy < 1) while the second is applicable when the deformations are dominated 
by the shear terms, i.e e xx <C and e yy <C ^ JL . The scheme without corrective term 



d B j _^^\d ( A j-l : B j\ _ B j =Q (16) 

txXi-Dxi n \i [f-^Vi < D Vi) D m U ' \ LV ) 



A convergent scheme for initial shape identification 6 

\t -jfl (A^r 1 + Hi) - e yy d (X/, 1 + B v {) - B y { + h.o.t. , (15) 

where (AJ-, Ay?) and (B^ , By?) are the x-y components of the residual vector A\ and 
the corrective term B\ respectively. If the corrective term vanishes, it is easy to show 
(see [8J) that two extreme cases should produce convergent schemes. The first one 

dy 

______ . , 1 I txrnilo f no copnnrl ic annlirio nlci «rnon f no r\ omrm oil 

n t* form o l a *t _______ , _____ 

dy c yy ^ 9x 

will be referred to as scheme (I) in the following. A suitable choice for the corrective 
term B\ should allow an unconditional convergence. One immediately sees that if B\ is 
chosen such that, 

dUx 
dy 
dTT 

■ -Jjf '/ (Xj ' + B4) - e yy fB y { -B y { = , (17) 

equations (fTlj) and (fTo]) simply reduce to 

A x j = -e^A^ 1 and A^ = -e w f A.f 1 , (18) 

which clearly defines a convergent scheme providing the absolute value of the strains is 
smaller than one which is necessarily the case in small deformations. Note however that 
since the convergence analysis is based on the Taylor expansion of the displacement field 
around Mf and that only the first order terms are retained, its validity is resctricted to 
situations when M{ is sufficiently close to Mf. Equations (fTEJ) and (JT7J) are rewritten 
in a more readable form as follows, 

(e d + l) B j + — \ d B 3 = - —\ d A j - 1 (19) 
\t XXl ^rijn Xl ^r \ t n yi , ^lyj 

^L\ d B j + (e d + l) B j - - ^v\ d A j ~ l (2Q) 
dx \ * ' * ~~ dx ' 

Since e xx f, e ra ^, ^-\f and ^r|f must be evaluated at Mf, these can be computed in 

practice at the first iteration and therefore the system of equations (JTUJ) and (|2U|) can 

be solved for the two components of the corrective vector B\. We denote by scheme 

(II) the one with the corrective term which satisfies equations (|19jl and (|2()jl exactly. 

Finally, a simple analysis of the orders of magnitude reveals that providing e xx d <C 1, 

e yy d ^ 1' T^flf ^ 1 anc ^ l^Hi ^ 1' ^he corrective vector reduces to B\ ~ 

A y ?~ l , — ^^-\fAj~ l ) and the latter scheme will be referred to as scheme (III). 

4. Test problem 

As a test case, a known and prescribed displacement field is considered and we seek the 
initial geometry so that the deformed one matches a disc of radius 0.01. The following 
displacement field is imposed: 

U = (U x , U y ) = (a(x + y), a(x - y)) . (21) 
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Figure 3. (a), (b) and (c): Initial and deformed geometries using the schemes 
(I), (II) and (III) respectively, a=0.6 and 100 boundary points. The arrows 
indicate the direction of increasing j. The figures also show the desired geometry 
and the true required initial geometry (analytical solution); (d): corresponding plot 
of max (| \Al 1 1) , i e [1,L] against the iteration number. 



Of course, in this case an analytical solution to the inverse problem can be found 
analytically by solving the following system of equations, 

x f - x ™ = a (x™ + yt ni ) , (22) 
yf - yT = a (xr - yf) , (23) 

for (x 1 ™, yl ni ) with (xf , yf) belonging to the circle of radius 0.01. This displacement field 
is however a good candidate to assess the proposed method since the terms e xx , e yy , 
and will all have the same magnitude equal to a and this is precisely the situation 
when, according to the previous analysis, the convergence of scheme (I) can not be 
guaranteed. 

Results are presented on Figure 01 for a=0.6, 100 boundary points [L = 100) and 
the three different schemes. Clearly, for the three schemes each iteration produces 
an initial geometry which comes closer and closer to the true required initial geometry 
calculated analytically using eqs. (j22J) and (j2HJ) and consequently the deformed geometry 
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Figure 4. (a): Initial and deformed geometries using scheme (II), a=0.9 and 100 
boundary points. The arrows indicate the direction of increasing j. The figure also 
shows the desired geometry and the true required initial geometry (analytical solution); 
(b): corresponding plot of max(||A^||) , i £ [1,L] against the iteration number. 



comes closer and closer to a disc of radius 0.01. The naked eye can readily distinguish 
that fewer iterations are necessary with scheme (II) (see Figure El (b)) to identify the 
required initial geometry than with the two other schemes (Figures El (a) and (c)). This 
is emphasized on Figure El (d) where the maximum of the norm of the residual vector 
from all the boundary points is plotted against the number of iterations. As predicted 
by the previous analysis, the convergence rate defined by the ratio of the norm of the 
residual vector at the current iteration to the one at the previous one, i.e. ||A|||/||A| -1 ||, 
is precisely equal to a=0.6 for scheme (II). This is simply a result of eqs. (fTHjl . It is 
worse for scheme (I) and scheme (III). The convergence rate is equal to 0.8486 in 
the former case and 0.84 in the latter. 

The true benefit of scheme (II) becomes even greater when the value of a is 
increased. As shown on Figure 0] (b), scheme (I) and scheme (III) fail to converge 
completely for a=0.9 whereas scheme (II) successfully converges towards the required 
initial shape. Again, in virtue of the previous analysis, the convergence rate is precisely 
equal to 0.9. The calculated initial and deformed geometries along with the desired one 
and the corresponding analytical solution are shown on Figure 0] (a). 

Finally, since scheme (III) is expected, according to the previous analysis, to be 
well suited for smaller deformations, an additional test is performed with a = 0.1. The 
results, plotted on Figure show a rapid convergence of the three schemes towards the 
required initial geometry. The convergence rate of scheme (III) is equal to 0.11 and is 
close to the value of 0.1 for scheme (II). Scheme (I) displays the worst convergence 
rate with a value of 0.1414 which confirms the benefit of adding at least an approximate 
corrective term in the updating of the boundary node location (eq. (J3J)). 
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Figure 5. Plot of max(||A^||) , i £ 
schemes and a=0.1. 



[1, L] against the iteration number for the three 



5. Conclusions 

This paper discusses the possibility to identify the required initial shape so that the 
deformed one matches very precisely a prescribed shape by means of an iterative scheme. 
It consists in using the desired work-piece geometry as an initial guess to the required 
initial geometry and back-tracking the latter by iteratively updating the locations of a set 
of boundary points. Based on an analysis of the convergence properties of the scheme, 
a corrective vector for the updating of the boundary point locations is derived. The 
addition of the latter is shown to extend the applicability of the scheme explored in jH] 
to arbitrary deformations when the shear and volumetric strains are equally important. 
The convergence is proven to be unconditional for small deformations and in the limit of 
the validity of the Taylor series expansion. Moreover, as illustrated by the test problem, 
significant improvement in the convergence rate is achieved by introducing this corrective 
term and to a lesser extent when the approximated corrective term is used. 

This scheme could offer a valuable alternative to other approaches for initial shape 
identification based on sensitivity analysis and optimization methods. At least four 
potential benefits may be outlined. No parameterization of the work-piece geometry 
is required since an arbitrarily large number of boundary points can be selected. The 
method is purely geometric and therefore its success does not depend on the type of 
constitutive law. The method can easily be used in combination with a commercial code 
for the computation of the displacement field, see [S] for an example. The question of 
how closely the deformed geometry matches the desired one which is not necessarily a 
trivial one is easily answered thanks to the introduction of the residual vector whose 
norm gives a clear measure of the mismatch between the geometries. This residual 
vector can only be defined because the set of reference points (at the desired locations) 
is defined at the first iteration of the scheme when the desired geometry is chosen as a 
first guess for the required initial geometry. 
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